Band touching from real space topology in frustrated hopping models 
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We study "frustrated" hopping models, in which at least one energy band, at the maximum 
or minimum of the spectrum, is dispersionless. The states of the flat band(s) can be represented 
in a basis which is fully localized, having support on a vanishing fraction of the system in the 
thermodynamic limit. In the majority of examples, a dispersive band touches the flat band(s) at a 
number of discrete points in momentum space. We demonstrate that this band touching is related 
to states which exhibit non-trivial topology in real space. Specifically, these states have support on 
one-dimensional loops which wind around the entire system (with periodic boundary conditions). 
A counting argument is given that determines, in each case, whether there is band touching or not, 
in precise correspondence to the result of straightforward diagonalization. When they are present, 
the topological structure protects the band touchings in the sense that they can only be removed 
by perturbations which aJso split the degeneracy of the flat band. 

PACS numbers: 71.10.Fd,71.20.-b,71.23.An,71.10.-w 
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I. INTRODUCTION 

The theory of "accidental" touching of energy bands 
in crystals has been recognized and studied since the 
early days of the quantum theory of solids. 1 By acci- 
dental, one means that the touching is not required by 
symmetry. A spectacular example of current interest is 
the Dirac point degeneracy of graphene, which leads to 
a host of interesting behavior^ Another class of heav- 
ily studied theoretical examples are the Dirac points ap- 
pearing in problems of two dimensional electrons mov- 
ing in periodic potentials in a magnetic field studied by 
Hofstadter 3 and others. Three dimensional Dirac points 
occur in models of unusual "spin Hall insulators" oc- 
curring with strong spin-orbit interactions A In all these 
cases, despite the accidental nature of the band touching, 
it is robust to perturbations of the Hamiltonian. This 
robustness has its origin in momentum space topology 
of the Bloch wavefunctions£ For instance, in graphene, 
each Dirac point is a source of a ±7r delta-function flux 
of Berry curvature, so that the line integral of the Berry 
connection, § c dk ■ Im.(u n fe|Vfe|u n fe) = ±7r for any curve 
C enclosing a Dirac point. If time reversal and inversion 
symmetries are maintained, the Berry curvature vanishes 
identically except at points of band crossing, and conser- 
vation of its flux protects the band crossings comprising 
the Dirac points. 

In this paper, we describe the topological protection 
behind a completely different instance of accidental band 
touching, which occurs in a broad class of "frustrated" 
hopping models. The models which we will consider actu- 
ally display in addition to band crossings a more dramatic 
phenomena: the presence of one or more completely Rat 
bands. Models with flat bands are particularly interest- 
ing physically because in this case the effect of interac- 
tions is wholly non-perturbative: interactions can recon- 
struct the states within the flat band manifold without 



any cost in kinetic energy. This is a powerful mecha- 
nism for generating complex and interesting many-body 
states, as attested by the richness of the fractional quan- 
tum Hall effect, which occurs as a result of the flat band 
degeneracy of Landau levels of electrons in a magnetic 
field. 

The frustrated hopping models we consider here arise 
in other contexts, e.g. the description of magnons in 
frustrated quantum antiferromagnets, and the motion of 
cold atoms in p-wave Bloch bands optical lattices. When 
the density of particles described by such frustrated hop- 
ping models is sufficiently low, short-range interactions 
of arbitrary (weak or strong) strength lead to particle lo- 
calization into a variety of crystalline patterns that are 
model specifici£' 7 -' & 9 ' 10 i 11 ' 12 These states are analogous 
to the Wigner crystal states of electrons in the Lowest 
Landau Level (LLL), which indeed occur for sufficiently 
small filling factor. These states can be understood as 
follows. From the flat band, one can construct single 
particle Wannier states (superpositions of wavefunctions 
with all momenta) which are strictly localized, i.e. have 
support on only a small finite number of sites. While 
Wannier states may always be constructed, only for the 
case of a flat band do they remain one-particle eigen- 
states. When interactions also have finite range and are 
repulsive, multiple particles can be present in spatially 
separated localized states with zero interaction cost, and 
the corresponding many-body wavefunction remains a 
many-body eigenstate. If the flat band has the minimum 
kinetic energy, such a state minimizes simultaneously ki- 
netic and interaction energy, and is therefore a ground 
state. Such states are generically possible for particle 
densities below some "close packing" threshold, at which 
interaction cost becomes inevitable. The states precisely 
at this threshold density are usually (but not always) 
periodic crystalline configurations. 

In the case of electrons in the LLL, the more inter- 
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esting fractional quantum Hall states occur above this 
filling factor, when there is some unavoidable interaction 
cost. However, because of the presence of a (large) gap 
between the LLL and the first LL, the interactions act en- 
tirely within the former. The zoo of fractional quantum 
Hall states is understood primarily from studies of the 
Coulomb interaction projected into the LLL subspace. 
A very interesting question is whether any similar rich- 
ness of behavior might occur, above the close packing 
threshold, when interactions are included in the frus- 
trated hopping models discussed here. This would appear 
a promising place to search for exotic orderless quantum 
spin liquid phases, loosely analogous to fractional quan- 
tum Hall liquids, that have been hypothesized to occur 
e.g. in frustrated magnet o 13 ' 14 . 

At this point the band touchings re-enter the picture as 
a hindrance to this search. The projection of the Hamil- 
tonian into the lowest energy flat band is strictly con- 
trolled only when there is a gap between this band and 
the higher dispersive ones, and when this gap is large 
compared to the strength of interactions. In the vast 
majority of frustrated hopping models (we will catalog 
many below), the gap actually vanishes due to touch- 
ings of the first excited band with the flat one at spe- 
cific points in momentum space. Before attempting to 
surmount this obstacle, it is crucial to know in each case 
whether this touching can removed by some small change 
in the Hamiltonian, or whether it somehow enjoys protec- 
tion that makes a search for such perturbations fruitless. 
The result of this paper is that in many cases the crossing 
is protected, and can only be removed by perturbations 
that also destroy the flatness of the low energy band. Like 
the protection of the Dirac points of graphene and others 
discussed above, the mechanism for this stability is topo- 
logical. However, because of the localized character of 
the states in the flat band, the topological structure lies 
in real space rather than momentum space. Specifically, 
the band touchings can be associated with eigenstates 
whose support is extended along non-contractible loops 
crossing a (toroidal) sample with periodic boundary con- 
ditions. 

The remainder of this paper is organized as follows. 
In Sec. HU we describe in detail the structure of local 
and topological loop states for one of our simplest exam- 
ples, the nearest-neighbor hopping model on a kagome 
lattice. We show how counting of these states requires 
band touching. In Sec. IIII1 we give a more abbreviated 
presentation of the generalization of these arguments to 
various other frustrated lattices. Finally, we conclude 
with a discussion in Sec. |lVl An appendix © describes 
an additional model with further-neighbor interactions 
without band touchings. 



II. KAGOME LATTICE MODEL 

The simplest model we will consider is the nearest 
neighbor tight hopping Hamiltonian on the kagome lat- 




FIG. 1: (Color online) Conventions for the shortest length 
Bravais lattice vectors for the kagome lattice (brown arrows) . 

tice, 

n t = -tY J (c\c ] + h. c ) , (i) 

m 

where the indices i,j denote the sites of the kagome lat- 
tice, (ij) denotes nearest neighbor pairs of sites, and the 
particles can be either Fermions or Bosons. We use this 
simple model to demonstrate all the generic flat band 
features discussed above. Aside of providing concrete ex- 
amples for all the unique features of flat bands, the anal- 
ysis of the kagome model also proves a good model with 
which to develop the techniques we will use throughout 
this manuscript. 

A. Band structure 

The band structure of Eq. consists of a single flat 
band with energy eo(q) = 2i and two dispersive bands 
with 

e±(q) = -<(l±V / 3 + 2A(q)) (2) 

where A(q) = cos(q • ai) + cos(q ■ a 2 ) + cos(q • a 3 ). Here 
ai,2,3 are the three shortest Bravais lattice vectors for 
the kagome (and triangular) lattice. Our conventions are 
described succinctly in Fig. [TJ The upper dispersive band 
e_(q) touches the flat band at the T point q = 0. 
The hopping term in the band basis is of the form 

n * = E / 4(qK(qMq) , (3) 

where the momentum integration is over the first Bril- 
louin zone, and we use v = 0, ± for the band index, and 
fi = 1,2,3 for the basis index. The new operators a„(q) 
are related to the original operators by a Unitary trans- 
formation. In particular, momentum eigenstates of the 
flat band consist of 

3 

ao (q) = 5>;(qMq) (4) 
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with tpfj,(q) 



sin(q-a M+2 /2) 



, where the greek index arith- 



ELi sin 2 (q ■ 



V(3-A(q))/2 

metic is always modulo 3, and 3 ~^( q ) 
a^+2/2). 

B. Localized states 



We can construct localized eigenstates by taking the 
linear combinations 



Al=Af 



5 -iq-R t 



4(q)V(3-A(q))/2, (5) 



with Af being some normalization. Here and elsewhere 
we will use A R to denote the creation operator for the 
localized eigenstates. Choosing R to be the position at 
the center of an hexagonal plaquette of the lattice, and 
normalizing the operator we find 



At 



R 



V .7 = 1 



Yd 



(6) 



where the indices 1 ... 6 enumerate the 6 successive sites 
around the hexagonal plaquette, as illustrated in Fig. [2] 
These local operators are very useful, but they are un- 
fortunately not canonical bosons or fermions. Rather, if 
Cj are bosonic, the commutation relations are 



^R'^R' 



= 5 



R,R' 



6 



(7) 



where the matrix Trr/ is the adjacency matrix of the 
triangular lattice formed by the centers of the plaquettes. 
For fermions, Eq. ([7]) holds with the commutator replaced 
by an anticommutator. 

The localized model can be understood directly in real 
space by considering a single triangle around the bound- 
ary of the plaquette. One of the corners has an amplitude 
of a second has ^= and a third has amplitude. The 
hopping amplitude from the first and second sites onto 
the third site cancels out. Thus the eigenstate is localized 
as a result of destructive interference, which is a very use- 
ful guiding principle in identifying these states in other 
flat band models. For a strictly localized wavefunction 
to be an eigenstate, the sum of hopping amplitudes onto 
sites outside the support of the wavefunction must vanish 
(see for illustration Fig.S]). 

One can create similar exact single-particle eigenstates 
on larger loops, by summing over the plaquette states on 
a number of contiguous plaquettes, and normalizing the 
state by the length of the boundary of the area covered 
by the plaquettes 



A ] - 



E 

ReA 



V6 



(8) 



Here A denotes the area covered by the plaquettes, and 
\dA\ denotes the length of the boundary of this area. In 
Fig. [5] we show one example of a three-plaquette loop. 




FIG. 2: (Color online) Depiction of localized eigenstates, on 
the boundary of a single and triple plaquette. Those sites 
with nonzero weight are denoted by a full (red) circle. The 
magnitude of the weights is always the same, but the phases 
alternate between ±1. The phases are denoted by ± signs 
next to the relevant lattice sites. 



C. State counting and band touching 

We now turn to the main question addressed in this 
paper, of the origin of the band touching. We will show 
that the set of localized eigenstates contains too many 
states to fit into the flat band alone. Specifically, the 
dimension of the space of localized state with the energy 
of the flat band has a dimension which is 1 larger than 
that of the flat band. This requires a contribution from a 
state of another band, which, since it is continuous, must 
touch the flat band at one point. 

Because the difference in question involves only a finite 
number of basis states (here 1, but there may be more 
in other examples in the next section), it is necessary to 
consider a large but finite system to make this count- 
ing precise. It is advantageous to use periodic boundary 
conditions (with a finite integral number of unit cells in 
each of two directions), since in this case the Bloch states 
in Eq. ((4]) remain eigenstates (with discrete q) in the fi- 
nite system. We must count carefully the number of lin- 
early independent states with energy eo. The plaquette 
states created by Eq. ([6]), naively all seem linearly inde- 
pendent, since they occur on different plaquettes. With 
open boundary conditions, the sum over all the plaque- 
ttes in the lattice leads to a state of the form of Eq. ([5]) at 
the boundary of the system. For periodic boundary con- 
ditions (putting the lattice on a torus), however, this sum 



vanishes since there is no boundary A q=0 = ^ R A R = 0. 
So when considering the Hilbert space spanned by the 
plaquette states © we have only (N — 1) independent 
states, where N is the number of plaquettes (and unit 
cells) in the lattice. This accounts for all but one state 
of the flat band. 
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FIG. 3: (Color online) The two non-contractible loop states 
around the handles of the torus. One loop consists of the 
sites marked by full (blue) circles, and the other by the empty 
(red) circles. As the other eigenstates in the flat band, the 
wavefunction has an alternating ± phase on the sites along 
the loops. 



The missing state is accounted for by a non- 
contractible loop around the torus. By decorating such 
a loop with alternating plus/minus signs, as illustrated 
in Fig. [31 one again satisfies the conditions for destruc- 
tive interference of outgoing waves, and the associated 
wavefunction represent an exact eigenstate, with the flat 
band energy. This state cannot be expressed as a sum 
of plaquette operators, or it would be possible to con- 
tract the loop just as any sum of plaquette states is. We 
have therefore found the missing state! However, we have 
an embarrassment of riches - there is not one such non- 
contractible loop, but two. In total we have (N + 1) 
states, all with the same energy. From the band struc- 
ture we know the flat band contains precisely N states, 
and so the additional state must come from another band, 
and for this reason one of the dispersive bands touches 
the flat band at exactly one point. 

In fact, from the loop states we can construct the plane 
wave Bloch state which touches the flat band explicitly. 
By taking an equal weight linear superposition of the 
non-contractible loops translated in any direction other 
than that along which the loop runs, one obtains a state 
with the same configuration in any unit cell, which there- 
fore has the Bloch form with momentum q = 0. The 
double degeneracy of states with q = signifies that not 
only must one of the dispersing bands touch the flat band 
at a point, but that the point is at q = 0. 




FIG. 4: (Color online) The localized states are exact eigen- 
states due to destructive interference between the hopping 
amplitudes from sites with nonzero weight (filled circles) to 
sites outside the boundary (empty circle). The lattice sites 
with nonzero weight are contained in a finite area, within a 
boundary marked by the dashed line. 



III. LOCAL EIGENSTATES 
A. Pyrochlore lattice model 

Taking the nearest neighbor hopping model 
(JXJ) on the pyrochlore lattice (instead of the 
kagome lattice) has two degenerate flat bands 
at eo = 2t, and two dispersive bands e± = 
-2t (1 ± y/1 + COS f COS f + cos f cos f + cos f COS f ) , 
where we have used the conventions ai = 4(0,1,1), 
a 2 = i(0, 1, 1), and a 3 = |(0, 1, 1) for the (FCC) Bravais 
lattice vectors, and eo = |(1, 1, 1), ei = |(— 1,1,1), 
e 2 = |(1, —1, 1) and e 3 = |(1, 1,-1) for the pyrochlore 
basis. Both flat bands touch the upper dispersive band 
at q = 0. The same localized plaquette modes that 
appear in the kagome model, are exact eigenstates for 
this pyrochlore model as well. However, whereas the 
number of hexagonal plaquettes in the kagome lattice is 
equal the number of unit cells, in the pyrochlore lattice 
the number of plaquettes is 4 times that of the number 
of unit cells. With two flat bands containing only 2N 
states, clearly these are not all linearly independent. 

Consider a volume enclosed by 4 plaquettes (see Fig[5]). 
Placing plaquette states with equal weight, and appro- 
priate relative signs, on each one of these 4 faces gives a 
total of zero. There are 2N such cells in the pyrochlore 
lattice, and therefore 2iV such constraints. This reduces 
the number of independent states we can construct out 
of the plaquette states to 27V. We choose to keep all the 
plaquette states for the plaquettes perpendicular to two 
out of the 4 (111) directions of the pyrochlore lattice. 

Now if we consider any one of the kagome planes along 
the two directions we chose above, we have the same 
additional constraint as in the kagome lattice - putting 
a plaquette state on every plaquette in the plane, with 
periodic boundary conditions results in zero, giving us 
one additional constraint. Taking into account the cell- 
constraints from the previous paragraph, there are only 
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FIG. 5: (Color online) Pyrochlore volume enclosed by 4 pla- 
quettes. Two of the 4 plaquettes are highlighted by thick 
(red and green) lines. Each one of the 2 plaquettes supports 
a localized eigenstate. 



two such linearly independent planes, so we have two 
additional constraints, reducing the number of linearly 
independent states we can construct from the plaquette 
states to 2N — 2. As in the kagome lattice, we now have 
the non-contractible loop states to consider. In the py- 
rochlore lattice there will be 3 such non-contractible loops 
rather than 2. Therefore, in total we have 2N + 1 states 
with the same energy of the flat bands. Exactly as argued 
in the kagome case, these particular non-contractible loop 
states can be made into q = states, and we therefore 
have 3 degenerate states at q = 0. The only way for the 
band structure to comply with this is to have one of the 
dispersive bands touch the two flat bands at the q = 
point which indeed is the case. 



B. Dice lattice model 

The nearest neighbor hopping model can give rise to 
a flat band on other lattices as well. We considered one 
example of a 2D lattice, and one example of a 3D lattice. 
In this subsection we shall mention one additional 2D 
lattice - the dice lattice, for a number of reasons. First, 
as opposed to the kagome and pyrochlore lattices, in this 
model the flat band touches dispersive bands at momenta 
other than q = 0. Second, it will be useful to compare 
between two different hopping models on this lattice, that 
both produce a flat band. The analysis essentially follows 
the same steps as in the kagome lattice model, and so we 
will not elaborate how the results were obtained. 

The dice lattice has a basis of 3 sites, two of which have 
coordination number 3, and one with coordination num- 
ber 6. In what follows, we shall refer to the latter sites as 
the coordination-6 sites. On the dice lattice, the nearest 
neighbor hopping model has one flat band at e = touch- 
ing two dispersive bands e± = ±2i\/2\/3 + 2A(q) (with 
A(q) the same as defined in the introduction) at the two 
momentum points q = ±(^=, 0) (with the Bravais lattice 



vectors taken with length 1). Adding an on-site potential 
which does not break the symmetries of the lattice(for in- 
stance an energy cost V to be on a 6-coordination site), 
one can gap one of the two dispersive bands away from 
the flat band, and only two degenerate points will remain. 
Therefore, our counting arguments will have accounted 
for N + 2 states with the energy of the flat band. 

The localized eigenstates of the dice lattice model are 
different than those of the kagome and pyrochlore lat- 
tices. Rather than residing in a loop around one or a 
number of plaquettes, the simplest localized states here 
have nonzero weight on the 6 sites neighboring a central 
coordination-6 site with alternating signs, as illustrated 
in Fig. El 

As for the kagome and pyrochlore lattices, a sum over 
all the localized states surrounding every 6-coordinated 
site, can produce zero with periodic boundary conditions 



A 



,,-±(4il o) = Er e+lq = 0. Apart from these two 
constraints, the localized states are all independent. The 
number of coordination-6 sites on the lattice is the same 
as the number of unit cells, and so we have accounted for 
N— 2 states. As in the kagome and pyrochlore models, we 
will find eigenstates composed of non-contractible loops 
around the torus. 

In Fig. [6] we show one of the two non-trivial loop states 
which exist for this model, with the weights being inte- 
ger powers of the factor ui = e ±l ~^~ , for a total of 4 non- 
contractible loop states. These are exact eigenstates pro- 
vided that 1 + uj + L)' 2 = is satisfied. Indeed, w — e 
are the two solutions of this equation. 



C. Honeycomb lattice p-band model 

Another hopping model with flat bands is the p-band 
hopping model on the honeycomb lattice introduced in 
Ref. i. In this model, only the planar p xy orbital 
states are considered at each lattice site. It is con- 
venient to describe any superposition of the two or- 
bital states on a site with an orbital unit vector vec- 
tor p = (jpx)Py) representing the state \p) = cos^|pa;) + 
sin0|pj,) with <fi the angle between the arrow and the 
x-axis. The hopping is assumed to occur only be- 
tween orbital states with the orbital vector parallel to 
the link. The resultant tight binding Hamiltonian is 



W = t E(y) Pi ■ ('< - ij) Pi ■ (r< - rj) 



-h.c. where the 

<1 



nearest- neighbor distance is taken as unity, and p^ ( Pj) 
are the creation (annihilation) operators for a particle at 
site j in the orbital state \p) . 

The two particle-hole symmetric flat bands in this 
model touch dispersive bands at q = 0. For the low- 
energy flat band, the local eigenstates are illustrated in 
Fig. [7] where the arrows denote the orbital state vectors. 
The high-energy flat band states have all the orbital- 
arrows pointing in the opposite direction on one sublat- 
tice, p —* —p. 

Exactly as for the kagome and pyrochlore lattice mod- 
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FIG. 6: (Color online) Dice lattice localized eigenstate de- 
noted by thick (red) closed loop, surrounding the site labeled 
by R. One of two non-trivial loops is indicated by a thick 
straight (blue) line. All sites with some particle weight on 
them are indicated by a filled circle. The amplitude is indi- 
cated on every site with non-zero weight. From this picture 
we can understand why the localized states are eigenstates, 
again invoking the picture of destructive interference. Con- 
sider the sites marked 1 and 2. Hopping from these sites can 
occur to either R or the site right outside the boundary of the 
localized state, in both cases, the hopping amplitude from 
sites 1 and 2 cancels out. Similar considerations for the other 
sites of the localized states yield the same result. For the 
non-contractible loop state, every site neighboring the sites 
with nonzero weight have 3 hopping amplitudes contributing 
1 + u) + lo 2 . As long as this sum vanishes, this is an exact 
eigenstate. 




FIG. 7: (Color online) p-band honeycomb local eigenstate and 
non-contractible loop state. The orbital states are denoted 
by a vector (in the present case the vector coordinates can be 
chosen real). The orbital vectors are always perpendicular to 
one link emanating from the site, and with the special form 
of hopping assumed in the model, the hopping amplitude on 
this link vanishes in this state. 



els, summing the localized states on all plaquettes results 
in zero &q =0 = X/R a R = 0- The localized states pro- 
duce N — 1 linearly independent states, leaving one state 
in each flat band unaccounted for. Two non-contractible 
loop states exist (one on each handle of the torus), with 
the same energy as the flat band states, their details il- 
lustrated in Fig. [7J A superposition of all the translations 
of a non-contractible loop state, results in a q = state. 
The two additional states then account for two q = 
modes in the band structure, and explain why the disper- 
sive bands must touch the flat bands at this momentum 
point. 



D. Dice lattice p-band model 

Another model that can be analyzed in a similar man- 
ner is the same p-band hopping introduced in Ref. 0, on 
the dice lattice. This model is different from the oth- 
ers presented here in that the flat bands are separated 
from all the dispersive bands by finite energy gaps. We 
will use the same framework to understand why the flat 
bands are gapped in this case. 

The band structure for this hopping model has two 
degenerate flat bands at e = 0, and two pairs of particle- 
hole symmetric bands with energies 

e(q) = ±tV2\J (6 + A)± ^/2(2A 2 -3A- 3B) , (9) 
with 

3 

A = J^cos (q- a M ) 

T ( l0 ) 

b = cos (q • ( a M + 2a P+i)) > 

where the three vectors ai,2,3 are the same minimal 
length Bravais lattice vectors indicated in Fig. Q] (now 
with length a/3) j making A, B two functions that are in- 
variant under the full symmetry group of the model. The 
dispersive bands are separated from the flat bands by a 
gap of A = \pk. 

Two local eigenstate modes can be found, on the same 
area unit surrounding a coordination-6 site. They are il- 
lustrated in Fig. [8] (type I) and Fig. [9] (type II), using the 
same conventions we have introduced for the honeycomb 
p-band model in the previous subsection. Note that the 
arrows indicating the orbital state are always in one of 6 
discrete directions, with the angles 0°, ±60°, ±120°, 180° 
from the x-axis direction. In these states, the arrow di- 
rections are always perpendicular to one link emanating 
from the site. 

Naively, all these states (type I and II) on different area 
units are linearly independent, resulting in 2N states that 
exhaust the number of states in the flat bands. However, 
as with the previous models explored in this paper, we 
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FIG. 8: (Color online) Type I p-band dice lattice local eigen- 
state. 




FIG. 9: (Color online) Type II p-band dice lattice local eigen- 
state. 

will find a number of constraints that show that this is 
not the case. 

Summing over the type I states on a set of non- 
overlapping area units that cover the entire plane results 
in zero, when we have periodic boundary conditions. An 
appropriate choice of the area units to sum over is il- 
lustrated in Fig. [lOl However, there are 4 such distinct 
sets covering the entire plane, and so there are 4 differ- 
ent sums (involving different sets of states) giving zero. 
The 4 sets of type I states are related by the Bravais 
translations of the lattice. 

As in the other models we have discussed, non- 
contractible loop eigenstates can be constructed, as il- 
lustrated in Fig. [TIJ The non-contractible loop eigen- 
states always come in pairs - one around each non-trivial 
loop on the torus. However, there are 4 pairs of distinct 




FIG. 10: (Color online) Non-overlapping set of area units. 
Summing over type I states on these area units adds up to 
zero. 




FIG. 11: (Color online) Non-contractible loop states. Only 
those sites with arrows on them have non-zero occupation. 



such states related by translation. The quadrupling is 
closely related to the 4 different unit area sets covering 
the plane discussed in the previous paragraph. The non- 
contractible loop states we are presenting here always lie 
on the edges of one of the 4 sets of area units. 

There are 4 additional constraints involving sums over 
type II localized eigenstates, and the non-contractible 
loop states. Adding the type II states on every area unit 
on an infinite length strip, comprising a non-contractible 
loop around the torus, we find the resultant state de- 
picted in Fig. [12J The resultant state is a sum of 4 non- 
contractible loop states, of the kind depicted in Fig. QTJ 
and so the type II states and the non-contractible loops 
are not linearly independent. The sum of type II on a 
strip translated by one Bravais lattice vector, perpen- 
dicular to the strip gives another independent relation 
between the type II states and the non-contractible loop 
states. Finally, for strips around the other handle on the 
torus, we find two additional such relations. 

In total, we have 27V localized eigenstates, 8 non- 
contractible loop states, and 4 + 4 = 8 vanishing lin- 
ear combinations of these 2N + 8 states. We therefore 
have precisely 2N linearly independent states, all with 
the same energy of the flat bands, exhausting the num- 
ber of states in the flat bands. Since no extra states at 
this energy need to be accounted for by the dispersive 
bands, a gap can occur in principle, and indeed shows up 
in practice. 
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FIG. 12: (Color online) Sum of type II states (marked by 
filled circles at their centers) on a non-contractible strip of 
area units. Only those sites with arrows on them have non- 
zero occupation. 



IV. DISCUSSION 
A. summary 

We have found that the presence or absence of band 
touchings can be understood for a variety of frustrated 
hopping models by a careful counting of linearly indepen- 
dent localized states. Crucial to this counting is the pres- 
ence, which we found for all the models examined here, 
of non-contractible loop states. In all cases, a vanishing 
superposition of the localized states was found, always 
involving a sum over area patches that cover the entire 
lattice, but some with different weights in the summation 
( non-zero q for instance). When working in a toroidal 
geometry, for which momentum q remains a good quan- 
tum number, the missing states eliminated by this van- 
ishing superposition were recovered as non-contractible 
loop states. The counting, for the case of a single flat 
band, is as follows. The total number of independent 
states with the flat band energy is N — M + dL, where 
N is the number of unit cells, M is the number of local- 
ized state independent superpositions that vanish (the 
number of missing states), and L is the number of dif- 
ferent 'flavors' of non-contractible loop states for each of 
the d handles of the torus (in d dimensions). When this 
is larger than N, the flat band is degenerate with some 
other band, at a finite number of points in momentum 
space. 



B. Are frustrated hopping models frustrated? 

We have referred to the hopping models discussed in 
this paper as "frustrated", simply because (apart from 
the p-orbital honeycomb example), they reside on lat- 
tices exhibiting strong geometrical frustration for anti- 
ferromagnetism. It is interesting to see if this abuse of 
terminology has any truth to it. Consider a general tight 



binding Hamiltonian with a flat band 

* ( n ) 

ij 

where the indices i,j include all the generalized coor- 
dinates of the particles (position, orbital state etc.). 
Choosing the particles to be fermionic, adding a spin | 
index and an on-site interaction term, we have 

H = Y1 c L*« c i« + U J2 c iT c iT c Jl c J'-l ■ ( 12 ) 

ija j 

Note that the interaction corresponds to the physical on- 
site repulsion only if the indices j represent only spatial 
coordinates, but this is only a minor complication we will 
ignore in this short and simplistic analysis. 

At half filling of the fermions, for infinite repulsive U, 
we get a Mott insulating phase, with one fermion oc- 
cupying each state of the (weighted) network of sites j 
described by the matrix tij. For finite but very large U, 
the virtual hopping of particles results in a Heisenberg 
model splitting the energy of the insulating states. It is 
very easy to show that the general Heisenberg model that 
emerges is 

H = ^J <i S i .S i + 0^ , (13) 

ij ^ ' 

At 2 - . i 

with Jij = —ff-, where Sj are the spin ^ operators. The 
matrix structure that is responsible for the flat bands 
carries over to the exchange interaction matrix. In par- 
ticular, if the hopping matrix takes on only two values, 
and t then = jjtij and the matrices simply differ by 
a multiplicative factor. 

When tij is non-zero for only nearest neighbor sites i, j, 
the spin model "descending" from the hopping Hamilto- 
nian is indeed geometrically frustrated. It is interesting 
to consider, however, a possibly more direct connection 
between flat bands and frustration. Indeed, somewhat 
rough arguments suggest that, fairly generally, the clas- 
sical spin Hamiltonian with exchange — jj-tij has an 
extensive ground state degeneracy when the lowest en- 
ergy eigenstates of tij form a flat band. We discuss some 
additional conditions below. 

The connection between flat bands and ground state 
degeneracy is through the Luttinger-Tisza method for 
finding ground states of classical spin models. The idea is 
the following. We first trade the normalization constraint 
on the spins, \Sj\ = S, for the weaker condition 

J2\Sj\ 2 = N s S\ (14) 
j 

where N s is the number of spins. With this weaker con- 
straint, it is simple to minimize H in Eq. (| 13|) . This is 
guaranteed to give an energy which is at least as low as 
for the minima of Ti taking proper spin normalization 
into account. The general solution is an arbitrary Teal 
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linear combination of the minimum energy eigenstates 
of Jij, and hence of tij. Because the hopping matrix is 
real, the eigenstates may also be chosen real. When there 
are rif flat bands, there are a total of n/N such eigen- 
states with minimum energy, for a lattice containing N 
unit cells. The solution is 

Nn f 

8^ = 5^8,^(0), (15) 

a=l 

where <frj(a) is the a th eigenstate of tij, and the s a is 
an unknown vector of real coefficients for each a. The 
total number of variables that may be varied is then 3 
real numbers for each a, and hence 3Nrif real numbers. 
Now we can attempt to impose the necessary constraints 
to get a physical minimum for the classical spin model. 
These consist of one constraint per spin on the spin mag- 
nitude. If the lattice contains rib spins in its basis (i.e. 
sites in the unit cell), then this gives Nrib constraints. 
Subtracting the number of constraints from the number 
of variables gives a total of (3n/ — rib)N degrees of free- 
dom remaining for physical minima of the classical spin 
model. Thus, for rib < 3n/, we are led to expect an ex- 
tensive degeneracy (macroscopic entropy) of spin ground 
states. This counting is certainly crude, and since the 
normalization constraints are non-linear, not entirely rig- 
orous. For instance, if rib = 3n/, it is probably the case 
that extensive ground state entropy may or may not be 
present, depending upon other details of the model. This 
is born out, for instance, by the case of the nearest neigh- 
bor kagome model, for which there is indeed an extensive 
ground state entropy, while rib = 3, rif = 1. However, we 
believe the conclusion provides a reasonable qualitative 
guide, though the estimate of the entropy density is prob- 
ably unreliable. Thus a large U Hubbard model with a 
kinetic energy with flat minimum energy bands indeed, 
when rib < 3n/, exhibits macroscopic ground state en- 
tropy, the classic signature of frustration. 

C. For the future 

The discussion in this paper is only a prelude to the 
study of interacting bosons and fermions in flat band sys- 
tems. As discussed in the introduction, this is understood 
in generality only for low density, below the appropriate 
"close packing" threshold. In this case the ground states 
are completely localized and minimize simultaneously the 
kinetic and interaction energies. Above this close packing 
density, the problem becomes much more intricate, and it 
is reasonable to expect delocalized "liquid" ground states. 
How to attack the problem in the range of densities for 
which the particles can still be accommodated in the flat 
band(s) but above close packing is an interesting open 
problem. It is intriguing to speculate that liquid states 
in this regime, at least for weak interactions, may have 
unconventional properties. With the increasing accessi- 
bility of such Hamiltonians in ultra-cold atomic systems 




FIG. 13: (Color online) Tight binding model with equal hop- 
ping amplitudes on three different link types. 

in optical lattices, clarification of this regime may well 
come experimentally rather than theoretically. 
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APPENDIX A: ANOTHER KAGOME MODEL 

Having identified a connection between frustrated 
magnetic models and flat bands, we examine one addi- 
tional model - one inspired by Ref. Hoi and mirro ring its 
geometric structure. The magnetic model in Ref. [l5j has 
been shown to support a spin liquid ground state. The 
magnetic model Hamiltonian consists of exchange inter- 
actions of equal strength on 3 different link types. This 
is described by an exchange matrix. We take this same 
matrix structure and construct a simple tight binding 
model with it. The Hamiltonian is precisely of the form 
of Eq. (fTTj) with the indices corresponding to the lattice 
sites of the kagome lattice. The hopping amplitudes are 
illustrated in Fig. [T31 and we will refer to this hopping 
model as the kagome-3 model. 

Our analysis finds the band structure of this model 
consists of two degenerate fiat bands at energy eo = +2t 
and one dispersive band ei = — t (4 + 2A(q)) with A(q) 
the same function introduced in the introduction. The 
single dispersive band has a minimum energy of — lOt 
at q = 0, and a maximum energy of — t at wavevector 
q = ^-x in our conventions. The flat bands are therefore 
gapped from the dispersive band with a gap of A = 3t. 
The convention we use for the 3 vectors illustrated in 
Fig. [U is ai = x, a 2 = — \x+ ^y, and a 3 = — \x — ^y. 

Proceeding as for the other models in this paper, the 
localized states we find are illustrated in Fig. [T5J and 
live on "bowtie" plaquettes. As in the pyrochlore lattice, 
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FIG. 14: (Color online) Illustration of the local constraint 
in the kagome-3 hopping model. Localized eigenstates living 
on the bowtie plaquettes are summed over the 6 plaquettes 
bordering a hexagonal plaquette in the kagome lattice. The 
figure shows two of these localized states, with the correct 
relative phase needed for the summation to vanish. Only 
those sites with non-zero weight are denoted by (red) circles 
(filled for +1, unfilled for -1). 




FIG. 15: (Color online) Localized eigenstates of the kagome- 
3 hopping model. The states have non-zero weight only on 
the 4 sites surrounding a 'bowtie' plaquette (the encircled 
regions). The sites with nonzero weight are marked by filled 
(red) circles, and with their relative signs indicated next to 
them. There are a number of different 'flavors' of bowtie 
plaquettes (and localized states) - 3 per unit cell. However 
there are only 2 independent bowtie states per unit cell. Here 
we show one choice of two independent bowtie states marked 
as A and B. 



naively there are 3 different "flavors" of bowtie plaquettes 
(and localized states). These can most easily be identified 
by considering how many bowties states involve a single 
up-pointing triangle (corresponding to a unit cell) . Were 
these states all linearly independent, we would have 3iV 
states in the flat band, rather than 27V (here as before, 
N is the number of unit cells). However, taking a bowtie 



plaquette state, rotating it around a hexagonal plaquette 
to produce 6 bowtie plaquette states around the hexagon 
( see fig. El), and finally summing these 6 wavefunctions, 
produces zero. For each hexagonal plaquette, of which 
there are N, there is one such constraint, and we find 
there are only 2 flavors of independent bowtie states - 
marked A, B in Fig. [TSl 

As in all the other examples we give in this manuscript, 
there are additional non-local constraints, which man- 
date the existence of non-contractible loop states. The 




FIG. 16: (Color online) Non-contractible loop states of the 
kagome-3 model. There are two types of the non-contractible 
loops. One kind are the exact same states appearing in the 
kagome model, and is denoted by a continuous (red) path 
between the sites with nonzero weight. The other kind of 
non-contractible loop states are denoted by dashed line, and 
consist of 2 sites with nonzero weight and alternating sign on 
each hexagonal plaquette in a chain of plaquettes. 



constraints we find involve vanishing summations over 
the bowtie plaquette states of type A ( see Fig. [15]), with 
wavevectors q = ib2,^(b!+b 2 ) (but not q = |bi). 
In the exact same manner, the sums over the type B 
plaquette states, with wavevectors q = ibi, \ (h\ + b 2 ) 
(but not q = ^b 2 ), also vanish. We therefore have 4 
constraints, and since there are a total of 2N flat band 
states, we expect 4 non-contractible loop states to exist. 
These can easily be found graphically, and are illustrated 
in Fig. HU 

As in section IIIIDl we find that in order for gapped 
flat bands to appear, we needed two different localized 
eigenstates occupying the same area unit, as well as non- 
contractible loop states that are not associated with one 
of the two sets of local eigenstates. 
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